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PREFACE 


The  research  described  here  was  performed  largely  under  Project  2303- 
F3-05  (Chemical  Structure/Bonding  Decomposition  Relationshipr)  and  Project 
2303-FA-03  (Molecular  Orbital  Calculations  of  Excited  Species).  A  portion 
of  the  funding  was  provided  by  the  AF  Armament  Laboratory  under  Project 


Order  ATL-1-68.  Mr.  Thomas  Floyd  was  the  AFATL  Program  Manager. 
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INTRODUCTION 


The  thermochemical  decomposition  of  2,4,6-trinitrotoluene  (TNT)  and 

related  nltroaromatlcs  in  the  liquid  phase  has  been  studied  by  a  variety 

of  techniques  for  many  years^  Janzen  first  observed  radical  formation 

during  the  thermochemical  decomposition  of  nitroaromatics  with  electron 

18 

paramagnetic  resonance  (EPR)  techniques  .  We  have  studied  the  thermal 

decomposition  kinetics  of  TNT  by  EPR^^.  This  paper  compares  EPR  studies 

19-22 

of  the  initial  TNT  radicals  with  MNDO  molecular  orbital  calculations 
on  TNT  and  model  systems  to  identify  the  initial  radicals  produced  and 
propose  plausible  first  steps  in  the  thermal  decomposition  of  TNT. 

EXPERIMENTAL 

EPR  spectra  were  recorded  on  a  Varian  E-109  spectrometer  operating  in 

26 

the  X-band  using  a  standard  TE^^q^  rectangular  cavity  as  previously  described 

A  custom-built  variable  temperature  apparatus  utilizing  a  Scanlon  Co.  heater 

maintained  temperatures  in  the  cafity  to  +  1°C  as  measured  by  a  Doric 
24  25 

trendicator.  TNT  and  trinitrobenzene  (TNB)  were  synthesized  using  pre¬ 
viously  described  procedures.  Hexamethylbenzene  (HMB)  was  purchased  from 
the  Aldrich  Chemical  Co.  All  compounds  were  purified  by  recrystallization 
from  ethanol  prior  to  use.  Mixtures  of  compounds  were  prepared  for  EPR 
analysis  by  mixing  the  components  in  the  molten  state  in  quartz  EPR  tubes 
at  temperatures  below  170°C. 

Simulation  of  EPR  spectra  for  comparison  with  the  experimental  spectra 
was  accomplished  with  a  program  written  and  generously  provided  by  Dr.  Ira 
Goldberg  of  Rockwell  Science  Center.  This  program  is  capable  of  simulating 
overlapping  spectra  from  multiple  radicals. 
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The  MNDO  and  MINDO/3  molecular  orbital  programs  were  used  to  calculate 
heats  of  formation  and  geometries  of  all  suspected  intermediates  and  tran¬ 


sition  states  for  mechanisms  tested.  A  related  force  constant  program, 
also  developed  by  M.  J.  S.  Dewar,  et.  al,  was  used  to  examine  suspected 
transition  states  and  provide  vibrational  frequencies  in  order  to  calculate 
tliermodynamlc  quantities.  Deuterium  kinetic  isotope  effects  for  some  steps 
also  were  calculated  using  these  thermodynamic  quantities. 

RESULTS  AND  DISCUSSION 

RADICAL  IDENTIFICATION: 

Previously  we  observed  at  least  two  radicals  from  the  thermochemical 
decomposition  of  pure  molten  TNT^^.  The  final,  or  product,  radical  has 
been  tentatively  identified  as  a  polymeric  radical  whose  number  of  radical 
site.s  grows  at  an  exponential  rate  throughout  the  decomposition.  The  initial 
radical,  whose  spectrum  is  composed  of  at  least  twenty  lines,  has  eluded 
identification  until  now.  By  careful  comparison  in  the  wings  of  the  spectrum 
of  pure  TNT  with  spectra  generated  by  doping  the  TNT  with  a  small  amount  o) 
tlMB,  we  find  that  the  TNT  spectrum  is  consistent  with  ttie  radical  structure 
below : 


Figure  1  shows  a  comparison  of  tlie  spectra  of  the  radical  formed  early 
in  the  decomposition  of  pure  TNT  at  d45°C  and  a  mixture  of  TNT  and  UMB 
(14:1  w:w)  at  ZOO^C.  The  spectrum  exhibited  by  the  HMB-doped  TNT  sample 


EARLY  TNT  RADICAL  SPECTRA 


Figure  1 


exhibits  a  large  coupling  constant  of  10.3  G  with  an  overall  intensity 

pattern  of  1:2:3:2:1,  indicative  of  two  equivalent  nitrogens.  Within  each 

of  these  five  groups,  a  six-line  pattern  is  evident  (coupling  constant  = 

2.7  G)  with  intensity  ratios  of  1:5:10:10:5:1,  indiiative  of  five  eqiix/a- 

lent  hydrogens.  The  computer  simulated  .spectrum  in  Figure  1  using  the 

measured  couplings  from  the  TNT/HMB  spectrum  shows  a  good  match  when  two 

I'quivalcnt  nitrogens  and  5  equivalent  hydrogens  are  assumed.  Ihe  spectrum 

matches  that  of  the  initial  radical  from  pure  TNT  decompos  i  t ion  in  Lite  wings. 

Tho  simplest  Interpretation  is  that  the  five  original  hvdrogens  present  i r. 

the  TNT  molecule  are  preserved  in  tlic  radical,  a.nd  one  nitre  group  has  been 

removed  leaving  two  equivalent  nitrogens.  The  apparent  equivalence  of  the 

ring  ami  metliyl  liydrogens  has  been  noti>d  for  other  radicals  fornied  from  TNT, 
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in  particular  tlio  radical  anion  .  It  is  also  interesting  to  note  that  the 
lIMI’i-doiied  sample  produces  the  initial  radical  spectra  at  much  lower  tempera¬ 
tures  tlian  pure  TNT  (-200'’C  instead  of  245°G)  and  the  .spectrum  remains 
practically  unchanged  for  over  thlrt>'  minutes  at  dOO^G. 

The  center  portion  of  the  neat  TNT  radical  spectrum  is  not  identical  with 
the  llMB-doped  TNT  spectrum,  even  though  the  wings  arc.  Simulations  which 
include  a  single  line  polvmer  radical  spectrum  with  the  llMB-d('|5ed  TNT  spec¬ 
trum  can  be  made  identical  with  the  I'.eat  TNT  spectrum  by  varying  the 
relative  amounts  of  the  two  radicals  (iiolymcr  +  TNI'  radical),  as  showTi  in 
Figure  2.  The  neat  TNT  spectrum  is  thus  app;iri'ntlv  the  superjuisi t ion  of 
tlte  polymer  radical  spectrum  on  that  of  the  TNI  r.ullcal  identified  in  t'u 
HMB-doped  experiments. 

A  reaction  analogous  to  the  IlMB-doped  TNT  decomros i t ion  apparently 
occurs  during  the  thermal  decompos  i  t  i  oti  <1!  1  ,  3 , 5-t  r  i  n  i  t  robt-nzt'no  (TNB) 


doped  with  a  small  amount  of  HMB.  TNB  Is  thought  to  be  very  stable  at 
temperatures  at  which  TNT  decomposes  rapidly  (.25n“C).  However,  the  HMli- 
doped  TNB  (TNB:  HMB  82:78  w/w)  produces  a  strong  radical  signal  at  200°c: 
which  is  quite  stable  as  shown  in  Figure  3.  An  overali  tive-line  pattern 
(coupling  constant  =  10.7  G)  with  a  1:2:3:2:1  intensity  pattern  again  indi¬ 
cates  two  equivalent  nitrogens.  Within  each  of  the  five  groups  is  a  four- 
line  pattern  (coupling  constant  =  2.6  G)  with  intensity  ratios  of  1:3:3:1, 
indicative  of  three  equivalent  hydrogens.  A  simulation  based  on  these 
values  shows  good  agreement.  Thus  TNB  doped  witli  HMB  appears  to  forrr.  a 
radical  analogous  to  the  TNT  radical  with  HMB: 


MNDO  CALCUI-ATIONS : 

Molecular  orbital  calculations  using  the  MNDO  and  MlNDO/3  methods 
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developed  by  Dewar,  ct^al,  have  been  carried  out  in  this  laboratory 

for  the  last  several  years.  Ti’.ese  methods  have  been  shown  to  be  useful 
in  a  variety  of  systems  in  the  determination  of  plausible  reaction  mecli- 
anisms  and  kinetic  parameters.  We  have  applied  these  metlieds  to  the  study 
of  the  initial  steps  of  TNT  and  related  nitroai  latics  thermochemica I 
decomposition.  Calculated  results  for  tiie  activation  energies  and 
deuterium  kinetic  isotope  effects  of  suspected  mechanistic  steps  are  com¬ 
pared  with  these  experimentally-determined  quantities  using  IDSC  and  EPR 


techniques.  Also,  structures  of  possible  Intermediates  are  compared  with 
experimental  spectra. 

Because  of  the  large  amounts  of  computer  time  to  completely  optimize 
the  geometries  of  the  TNT  species  involved  in  suspected  decomposition  steps, 
we  cliose  to  use  1-nitropropene  as  a  model  system: 


TNT  1 -n i t ropropene 


NO2 


Previous  experimental  evidence  had  indicated  involvement  of  a  methyl  group 
hydrogen  in  the  rate-determining  step.  Based  on  this,  six  potential  rate- 
determining  steps  were  selected  for  calculations  of  /LH'’  and  : 


fntramolccular  H  Transfer 


H 


Intramolecular  0  Insert  ion 


(i)  2  (TNT) 


Intermolecular  H  Transfer 
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ONO 


NO2 


Step  (5)  is  a  removal  of  an  NO2  group  (either  para  or  ortho).  Step  (6)  is  a 
nitro  group  rearrangement  to  form  a  nitrite  species. 

In  order  to  test  the  usefulness  of  the  1-nitropropene  model  in  simulat¬ 
ing  TN'I'  behavior,  AH°'s  for  all  six  steps  were  calculated  for  both  TNT  and 
the  1-nltropropene  model  system,  all  geometries  being  optimized.  The  results 
in  Table  I  are  MNDO  results.  MlNDO/3  results  gave  the  same  agreement  between 
the  two  molecular  systems  that  the  MNDO  results  did. 

TABJ.£  _I 

Results  of  MNDO  Calculations:  Reaction  AH°'s 


AH"  (TNT  System) 

AH°  (l-nitropropene  system) 

Stejp 

(Kcal/mole) 

(Kcal/mole) 

(1) 

I nt  ramolecular 

H 

Trans  f  er 

+  16.6 

+14.5 

(2) 

Intramolecular 

0 

Insertion 

-49.9 

-49.3 

(3) 

I ntermolecular 

H 

Transfer 

+  16.8 

+19.5 

(4) 

Intermolecular 

0 

Insertion 

-52.4 

-48,5 

(5) 

NO^  Removal 

+32.3 

r+26.9  (ortho-N02) 

(+32.5  (rara-N02) 

(6) 

Nitro  Rearrangement 

-38.4 

-41.1  (para-N02) 

Thus  It  appears  that  the  nitropropene  system  Is  at  least  a  potential 
model  for  simulating  TNT  decomposition  steps.  It  is  a  good  model  if  the 
calculated  activation  energies  arc  also  similar. 

We  simulated  the  reaction  pathway  for  each  of  tiie  six  steps  liy  defining 
a  reaction  coordinate,  fixing  it  at  various  intermediate  points  along  the 
pathway,  and  then  allowing  all  other  geometric  variables  to  optimi7e  at 
eacli  value  of  the  reaction  coordinate.  If  a  maximum  vas  located  along  the 
patiiway,  the  transition  state  optimization  program  was  run.  Final ]v,  tlu 
vihrational  force  constant  program  was  employed  on  the  transition  state 
g.eometry  to  support  the  existance  of  a  transition  state  by  the  presence  o: 
one  i)eg,.ativo  vibrational  force  const.ant,  which  is  idi'nlifiiil  as  the  reaction 
coord  In.ite .  Kach  of  the  six  steps  will  be  discussed  in  turn. 

The  intramolecular  hydrogen  transfer  (1)  calculations  .''or  1 -ii  i  i  ropropene 
produced  a  species  which  proved  to  be  n  transition  ^Male  witti  an  .activation 
energy  of  61.8  keal/mole.  Tlie  most  recently  measuiaal  valui'  (FPR)  for  the 
•  ictivatLon  is  40.9  kcal/mole^  ^ .  The  optimi.aed  geomi  t  rv  from  tlu-  calculations 
was  intermediate  between  that  of  the  reactant  and  product,  with  the  0-H  bond 
order  being  0.48  and  the  r-H  bond  order  fl.ld.  flu-  geora  t  lu  is  slii'wn  in 
Figure  4.  A  c;j  1  eu  1  a  t  i  on  of  a  deuterium  kinetic  i'-otirpe  I'Mi-ci  iitili/inc,  the 
ca  I  eu  la  t  ed  force  ennstant-  .iiul  miime:\ts  of  inertia  vieided  .i  k,,/k,  value  I'f 
d.d  at  227°C.  The  measured  value  of  k.j/k,,  iov  tlie  r.tt  -.‘e  f  orm  i  n  i  ng  st  t  p  ol 
TNT  docomp<’si  tion  near  this  same  t  emin-ra  t  ure  i 1.'*  liiu-  liie  call  illa¬ 
tions  do  not  rule  out  this  intramolecular  livdri  an  transfer  as  the  rate- 
determining  step  for  TNT  thermoehemic.il  deeompos  i  t  i  on  ,  particularli  sin.  i 
MNDO  usu.ally  overestimates  activatii'c.  energies  !et  hydrogen  transler  reae- 


t  ions . 


TRANSITION  STATE  FOR  INTRAMOLECULAR 


HYDROGEN  TRANSFER 


Figure  4 


TRANSITION  STATE  FOR  INTERMOLECUL AR 
OXYGEN  INSERTION 


ARE  HYDROGEN 

Figure  6 


On  the  other  hand,  the  intramolecular  oxygen  insertion  (2)  calculations 
for  ] -nitropropene  gave  very  high  energy  barriers  even  though  AH®  for  this 
step  is  quite  negative.  In  one  case  for  which  tlic  methyl  group  was  not 
allowed  to  rotate  while  the  oxygen  was  moved  toward  the  carbon,  the  molecule 
dissociated  into  acetylene  and  methyl  nitrite  witli  an  activation  enegy  of 
96  kcal/mole.  If  complete  relaxation  of  the  mctiiyl  group  was  allowed,  no 
stable  product  was  formed.  Thus  it  appears  that  this  step  is  not  likely 
during  the  thermochemical  decomposition  of  TNT. 

The  intermolecular  hydrogen  transfer  (3)  transition  state  has  not  yet 
been  located,  but  approximations  point  to  an  activation  energy  of  82  kcal/ 
mole.  This  step  also  appears  to  have  too  high  an  activation  energy  to  be 
a  viable  candidate  for  the  first  step  of  TXT  decomposition. 

'I'hc  intermolecular  oxygen  insertion  reaction  (4)  produceil  a  transition 
state  with  an  an  activation  energy  of  110  kcal/mole.  The  geometry  of  this 
transition  state  is  shown  in  Figure  5. 

Wo  also  investigated  rearrangement  of  an  NO,  group  and  removal  of  an 

U 

NO^  group  from  the  molecule,  even  though  these  steps  vtould  not  produce  a 

primary  deuterium  kinetic  isotope  effect  as  indicated  by  the  previous 

15 

experimental  evidence  .  We  found  that  tiie  calculated  dissociation  energy 
for  1-n  itropropene  of  32.3  kcal/mole  agreed  well  with  the  experimental 
value  for  the  activation  energy,  particvilarly  since  MXDO  typically  under¬ 
estimates  C-NO^  bond  energies.  In  addition,  Liic  calculated  deuterium 
kinetic  isotope  effect  at  227”C  for  NO.,  removal  was  1.7,  even  though  it  is 
a  secondary  kinetic  Isotope  effect.  This,  too,  agrees  well  with  the  experi¬ 
mental  value  of  approximately  1.7.  The  nitro  group  ’-carrangemont  to  nitrite 
for  l-nitropropene  had  a  calculated  activation  energy  of  52.4  kcal/mole. 


slightly  more  than  the  experimental  number  for  TNT.  The  geometry  of  the 
transition  state  is  shown  in  Figure  6.  This  step  cannot  be  ruled  out  as  a 
potential  first  step.  In  addition,  if  NO  then  dissociates  from  the  molecule, 
the  resulting  radical,  with  an  oxygen  on  the  para  position,  could  give 
qualitatively  the  experimentally  observed  EPR  spectrum. 

Table  IT  summarizes  activation  energies  and  deuterium  kinetic  isotope 
effects  for  both  1-nitropropene  and  TNT.  The  activation  energy  given  for 
TNT  for  tlie  Intramolecular  hydrogen  transfer  step  is  a  lower  bound,  and 
the  one  for  the  intermolecular  hydrogen  transfer  is  an  upper  bound. 

TABLE  II 

RESULTS  OF  MNDO  CALCULATIONS: 

ACTIVATION  ENERGIES  (KCAL/MOLE)  AND  DEUTERIWI  ISOTOPE  EFFECTS 


I’ROCESS 

E  (1-NP) 
a 

E  (TNT) 
a 

K^/K^(1-NP) 

(1) 

Intramolecular 

H  Transfer 

63.8 

>55 

2.2 

(2) 

Intramolecular 

0  Insertion 

96 

O) 

NO 2  Removal 

32.3 

^26.9(o-NO,) 

1.7 

(4) 

Nitro  Rearrangement 

52. 

_ _ 

_ 

(5) 

I n  t  ermo 1 cc  u la  r 

H  Transfer 

.82 

<122 

_ 

(6) 

1 ntermol ocular 

0  Insertion 

110 

CONCLUSIONS 


An  inlti.al  radical  in  the  thermochemical  decomposition  of  TNT  and 
related  nitroaromat ics  is  one  in  which  a  nitro  group  or  an  NO  group  has 
been  removed  from  the  TNT  molecule.  MO  calculations  for  removal  of  a 


nitro  group  from  the  molecule  are  consistent  with  experimental  evidence 


both  for  the  activation  energy  and  the  deuterium  kinetic  isotope  effect. 

The  enhancement  of  radical  production  upon  addition  of  H>iB  indicates  th.at 
the  interaction  of  a  second  species  results  in  much  easier  NO^  removal. 

'i'wo  likely  mechanisms  arise: 

(1)  The  donor  species  (HMB)  transfers  a  hydrogen  atom  to  the  acceptor 
species  (TNT  or  TNB) .  The  removal  of  HONO  can  then  be  possibly  accomplished 
more  easily  than  the  original  nitro  group. 

(2)  A  charge-transfer  complex  (ffMB/TNT  or  HMB/TNB)  is  formed  which 
weakens  in  some  manner  the  C-nltro  group  bond.  Thus  NO^  can  be  removed 
more  easily. 

The  calculations  indicate  that  oxygen  atom  insertions  are  unlikely  due 
to  the  high  energies  of  the  transition  states  involved.  Also,  the  activa¬ 
tion  energy  for  Intermolecular  hydrogen  transfer  has  a  value  which  is 
probably  inconsistent  with  the  experimental  value.  Intramolecular  hydrogen 
transfer  or  nitro  group  rearrangement  to  the  nitrite  cannot  be  ruled  out  as 
tlie  initial  step,  although  they  appear  to  be  less  likely  than  simple  nitro 
group  removal. 
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ABBREVIATIONS,  ACRONYMS.  AM>  S.YMBOl  S 


’I'NT  -  2,4,6  Trinitrotol'^enc 

F.PR  -  iMoctron  Paramagnetic  Resonanci- 

MNDO  -  Modified  Neglect  of  Diatomic  Overlap 

MINDO/1  -  Modified  Intermediate  NegK’ct  of  D  i  f  fet  i-nt  i  a  1  C'verlap/l 
-  Transverse  Electric  Mode  10? 

TNB  -  1,1,5  Tr initrohenzene 
lIMB  -  Ilexamethy  1  henzene 
C  -  Canss 

w:w  -  Weight  to  Wi-ight  Ratio 
°C  -  Degrees  (Centigrade 

AH'’  -  Change  of  enthalpy  for  (he  react  ion 

E.  -  Ai-tivation  energy  for  the  reaction 
a 

1DS(C  -  Isothermal  Differential  Scanning  Calor imct  i  i 
X  -  Nitro  group 

Ki|/Kd  -  Deuterium  Kinetic  isotope  ei'ei  t 


